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ABSTRACT 

Aims. Stellar activity is a potential important limitation to the detection of low mass extrasolar planets with indirect methods (RV, 
photometry, astrometry). In previous papers, using the Sun as a proxy, we investigated the impact of stellar activity (spots, plages, 
convection) on the detectability of an Earth-mass planet in the habitable zone (HZ) of solar-type stars with RV techniques. We extend 
here the detectability study to the case of astrometry. 

Methods. We used the sunspot and plages properties recorded over one solar cycle to infer the astrometric variations that a Sun-like 
star seen edge-on, 10 pc away, would exhibit, if covered by such spots/bright structures. We compare the signal to the one expected 
from the astrometric wobble (0.3 /jas) of such a star surrounded by a one Earth-mass planet in the HZ. We also briefly investigate 
higher levels of activity. 

Results. The activity-induced astrometric signal along the equatorial plane has an amplitude of typ. less than 0.2 ^as (rms=0.07 //as 
), smaller than the one expected from an Earth-mass planet at 1 AU. Hence, for this level of activity, the detectability is governed by 
the instrumental precision rather than the activity. We show that for instance a one Earth-mass planet at 1 AU would be detected with 
a monthly visit during less than 5 years and an instrumental precision of 0.8 fias . A level of activity 5 times higher would still allow 
such a detection with a precision of 0.35 yuas . We conclude that astrometry is an attractive approach to search for such planets around 
solar type stars with most levels of stellar activity. 

Key words. (Stars;) planetary systems - Stars: variable - Sun: activity - (Sun:) Sunspots - Astrometry. 



1. Introduction: 

Stellar activity is now recognized as a potentially strong limita- 
tion for the indirect detection of planets. Indeed, spots and bright 
structures (plages, network) produce brightness inhomogeneities 
at the stellar surface that affect the photometric, astrometric and 
radial velocity (RV) signals. The RV signal is also affected by 
the inhibition of convection in the active area. The amplitude of 
the activity-related stellar noise depends on the activity pattern 
and intensity, which is related to the stellar properties (age, tem- 
perature, etc). In the cases of young, active stars, or late-type 
stars, the signal may mimic that of a giant planet with periods 
of the order of the star rotational period. In the case of less ac- 
tive, solar-type Main-Sequence stars, the noise is much lower, 
but can nevertheless affect the detection of terrestrial planets. 
In two recent papers, we have investigated the impact of spots 
(Lagrange et al. 2010 herefater paper I), plages and convection 
(Meunier et al. 20101 herefater paper II) on the detectability of 
an Earth mass planet located in the habitable zone (HZ) of the 
Sun, as seen edge-on and observed in RV or in photometry. To do 
so, we took into account all spots and plages recorded over one 
solar cycle. We showed that providing a very tight and long tem- 
poral sampling (typically twice a week over more than 2 orbital 
periods) and an RV precision in the 10 cm/s range, the photo- 
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metric contribution of plages and spots should not prevent the 
RV detection of Eaith-mass planets in the HZ. This is no longer 
true when convection is taken into account. 

Here, we complete the previous work by estimating the as- 
trometric signal produced by the same structures over the same 
solar cycle, and we compare it to the astrometric signal induced 
by an Earth-mass planet located at 1 AU from the Sun. We also 
compare this signal to the RV and photometric signals. 

2. Data and results 

2.1. Data 

The data used to compute the astrometric signals are fully de- 
scribed in paper II. In brief, they consist of all (20873) sunspot 
groups larger than 10 ppm of the solar hemisphere, provided 
by USAF/NOAA (http://www.ngdc.noaa.gov/stp/SOLAR/) and 
all (1803344) bright structures : plages (in active regions) and 
network structures (bright magnetic structures outside active re- 
gions) found from MDI/SOHO ( iScherrer et al. 199 5) magne- 
tograms larger than 3 ppm, during the period between May 5, 
1996 (JD 2450209) and October 7, 2007 (JD 2454380). The tem- 
peratures associated to these two types of structures was deduced 
from the comparison between the variations of the resulting iiTa- 
diance and the observed ones (Fr oehlich and Lean 19981 1 during 
this period. Our temporal sampling is one day, with a few gaps, 
corresponding to a coverage of 86 %. 
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Table 1. Measured rms for the shifts, RV and photometry due to 
the spots and bright structures 



Period 


rms(AX) 


rms(AY) 


rms(RV) 
w/o conv. 


rms(RV) 
with conv. 


rms(TSI) 


all 


0.07 


0.05 


0.33 


2.4 


3.6 10-" 


highl 


0.09 


0.06 


0.42 


1.42 


4.5 10-" 


high2 


0.08 


0.05 


0.37 


1.62 


3.9 10-" 


low 


0.02 


0.01 


0.08 


0.44 


1.2 10-" 




1000 2000 3000 4000 
Time (JDB-2450000) 
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Notes. Shifts in x and y are in /jas , RV are in m/s. RV and photometry 
rms are taken from Paper 11. The entire cycle (refered as "all") as well 
as low and high activity periods are considered. The low activity pe- 
riod is from July 1, 1996 (JD 2450266) to April 1, 1997 (JD 2450540), 
and the first considered high activity period (refered as "highl") is from 
February 1, 2000 (JD 2451576) to November 1, 2000 (JD 2451850). 
Both correspond to the low and high activity 9-months periods consid- 
ered in Paper II. We also consider a longer (923 days) activity period 
(refered as "high2"), from JD 2451577 to JD 2452500. 



Fig. 1. Temporal variations of the astrometric shifts along (top) 
and perpendicular to (bottom) the equator, due to the combina- 
tion of spots and bright structures. The Sun is supposed to be 
seen edge-on and located 10 pc away. 
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Fig. 2. Astrometric shift (deltax) as a function of time (top) and 
excursions of the photocenter (bottom) during a low activity pe- 
riod (JD 2450335 to JD 2450365). 



2.2. Astrometric temporal variations 

For each date, we compute the position (x,y) of the Sun photo- 
center due to these dark and bright structures, assuming it is seen 
edge-on, 10 pc away. The resulting shifts along and perpendic- 
ular to the equatorial plane are shown as a function of time in 
Fig.[T] We report in Table [T]the rms of the shifts during the entire 
cycle, as well as during low activity and high activity period, and 
the corresponding values for the RV and photometric signals. 

Finally, we provide in Fig. |2] the temporal variations of the 
astrometric shift along the equator as well as the global excur- 
sion of the photocenter during a short 30-day period during the 
low activity phasis. The astrometric signal, if measured precisely 
enough, allows the determination of the star rotation period. So 



would also the RV and photomteric signals. The astrometric sig- 
nal would also provide additional and unique information. In 
Fig.|2j we see that the structure at the origin of the observed vari- 
ations is located on the southern hemisphere of the star, it mov- 
ing from East to West, ie from the left to right in the lower panel 
of Fig. |2] Therefore, the astrometric signal can allow, within 
one rotation period, to determine the direction of the star ro- 
tation axis in the plane of the sky. Note that this assumes that 
we know by photometry or RV measurements if the perturbat- 
ing structures are hotter or cooler than the photosphere and that 
this requires a very good astrometric precision. Also, one needs 
to isolate and monitor one single structure, so only periods of 
low activity can provide such information (indeed, during more 
active periods, the contributions of all individual structures mix 
up). 



2.3. Relations between RV, photometric and astrometric 
signals 

The astrometric variations along the equatorial plane are well 
correlated with the RV signal due to spots and plages (see top 
panel of Fig. |3]l. For a single spot, the ratio between the shift 
and the RV depends on the star's properties (v. sin /) and is 
almost independant from the spot properties. The correlation 
found here shows that this property remains when taking sev- 
eral spots/structures into account. The slope of the RV variations 
(in m/s) to the shifts (in juas ) is 5, compatible with the value 
expected from a single spot. Once convection is taken into ac- 
count in the RV computation, however, there is no correlation 
any longer between the astrometric shift in x and the RV varia- 
tions (see Fig.|3]l. This difference of behaviour between the "no- 
convection" and "convection" cases could in principle be used 
to estimate the level of convection, if any, associated to the star. 

Finally, no correlation is found between the relative pho- 
tometric variations and the astrometric shifts (either the Ax or 
the shift to center). This is not surprising as the photometric 
variations depend marginally on the spots/plages locations, con- 
versely to the astrometric shifts. 
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Fig. 3. Astrometric shifts along the equatorial plane versus RVs 
due to spots and plages (top) and to spots, plages and convection 
(bottom). No instrumental noise is assumed here. 



3. Discussion 

3. 1. Planet detection in astrometry 

We see in Fig. [T] that the activity-induced shifts are most of 
the time smaller than 0.2 yuas, ie smaller than the amplitude 
of the shifts induced by an Earth mass planet at 1 AU (0.33 
juas along the equatioral plane and //as perpendicular to the 
equator). These values are in agreement with the estimations of 
[Makarov et al. 20101 The rms of the equatorial shift over the en- 
tire cycle is ^ 5 times smaller than the amplitude of the Earth sig- 
nal, and 4 times smaller during high activity. This means that the 
Sun activity would not prevent to detect an Earth located in the 
HZ, provided the precision of the data allows such mesurements. 
This is illustrated in Fig. |4] where we show the periodogramme 
of the astrometric shifts along the equatorial plane once a one 
Earth-mass planet has been added (no noise). In this example, 
we use a limited amount of data instead of the whole set of data, 
in order to be closer to a real case. The chosen temporal sam- 
phng is of 1 month t 5 days over 50 months, typical of the strat- 
egy adopted for the NEAT instrument recently proposed to ESA 
(Leger2010i, dedicated to a systematic astrometric search for 
exo-Earths in the HZ of nearby stars. 

We then explored different instrumental noises. An exam- 
ple is given in the lower panel of Fig. |4] where we assume a 
noise of 0.8 juas The temporal sampling is the same as be- 
fore. The planet is still detectable. However, in that case, the 
coiTelation seen in the previous section between the astrometric 
shift along the equatorial plane and the RV variations is not any 
longer present. This is due to the fact that the instrumental noise 
is larger than the activity-induced noise. The consequence is that 
if we take the instrumental noise into account, for this level of 
activity, it will not be possible to estimate the level of convection 
using the astrometric and RV data. 

We finally explore different levels of activity. We found that 
an instrumental precision similar to that expected for the NEAT 
instrument for a 5 hours-long visit of a G2 star located 10 pc 
away (0.35 //as ) will easily allow the detection even for stars 



' this value corresponds to the precision that NEAT will provide in 1 
hour observation. In practice, each visit will be longer (up to 5 hours) 
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Fig. 4. Periodogrammes of the shift along the equatorial plane 
of the Sun viewed edge-on and surrounded by a one Earth-mass 
planets at lAU, and observed once a month approximately, dur- 
ing 50 months. From top to bottom: 1) low activity peiiod, no in- 
strumental noise; 2) high activity period, no instrumental noise; 
3) low activity period, 0.8 //as rms instrumental noise; 4) high 
activity period, 0.8 //as rms instrumental noise. 




Fig. 5. Periodogrammes of the shift along the equatorial plane 
assuming a Sun 5 times more active than observed, viewed edge- 
on and suiTounded by a one Earth-mass planets at 1 AU, and 
observed once a month approximately, during 50 months. A 0.35 
//as rms instrumental noise has been added to the data. Left: low 
activity period; right: high activity period. 

5 times more active than the Sun (we assume conservatively 
that such a star will produce shifts five times higher than the 
Sun) This is illustrated in Fig. [s] where we show the peri- 
odogrammes of the shift along the equatorial plane (same sam- 
pling as before) during the low and high activity periods. Given 
ILockwood et al. 2007 1 relation between the relative photometric 
variations and the Ca H & K line strength index log(R'HK), and 
given the log(R'HK) for the Sun (between ^ -4.85 and -5.0 
respectively during the high and low activity phases), this cor- 
responds to a log(R'HK) smaller than -4.25. Most of active 
stars fall into this regime. Indeed, among the 385 Hipparcos stars 
closer than 20 pc with a known R'HK, only 8 have log(R'HK) 
> -4.2 (see Fig.|6j Maldonado et al. 2010)). 

So far we have assumed that the angle of the planet orbit 
(supposedly identical to the equatorial plane) is known. In the 
case of real observations though, this information will not be 
known a priori and the planet may orbit outside the star's equa- 
torial plane. It appears that even in the presence of stellar activity 
and instrumental noise at a level of 0.35 //as per measurement, 
the astrometric data will allow to recover this information. We 
consider the example of a planet on an orbit inclined by -60 de- 
grees with respect to the North. The star is supposed to have ei- 
ther the same activity level as the Sun, or 5 times this level. The 
astrometric shifts are measured in a referential which is allowed 
to rotate from (equatorial plane) to 360 degrees. We measure 
in each case the shifts projected on this referential, and the asso- 
ciated rms along each axis of the referential. We consider either 



to ensure a better precision (precision is proportional to V'obs) 



- in practice, a level of activity five times higher will probably in- 
duce shifts smaller than 5 times the shifts, as statistically, there will be 
more chances that the effects of structures located on a given side of the 
hemisphere be balanced by that those of structures on the opposite side 
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Fig. 6. Histogram of log(R'HK) for all Hipparcos stars closer 
than 20 pc, with a known R'HK. 
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Fig. 7. Rms of the shift along the x-axis of a rotating referential 
as a function of its rotational angle (see text). In red, case of solar 
activity. Jn blue, case of 5 times the solar activity. In all cases, 
we consider that the instrumental noise is (plain lines) or 0.35 
microas (dashed lines). The vertical dashed line shows the angle 
of the inclination of the planet orbital plane. 



24, 4 and 14, respectively. Therefore, for RV techniques, the de- 
tection limit is set by the stellar activity (and mainly by con- 
vection effects) rather than the instrumental precision (provided 
precisions of 10 cm/s are available). 

To prepare and accompany a space astrometric misson ded- 
icated to Eaith-mass planet detection, a number of previous RV, 
or photometric monitoring would certainly be very precious. In 
particular, RV measurements are very powerful to derive the 
properties of the star's activity. They are in fact more power- 
full than astrometry as the rms of the activity-induced RV varia- 
tions is much larger than the instrumental precision (conversely 
to astrometry). Previous long-term RV monitoring can allow to 
identify the presence of larger bodies in the system and char- 
acterize them (note that lower precision astrometric monitoring 
can also do this). During the astrometric measurements, simulta- 
neous RV observations with a precision of about 10 cm/s would 
provide detailed information on the actual activity of the star and 
will help the analysis of the astrometric data. 
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no or 0.35 juas level instrumental noise, and either a solar- like 
level of activity or 5 times this level. The temporal sampling is 
the same as before (50 data points). We show in Fig.[7]the rms as 
a function of the rotation angle of the referential. We see a clear 
maximum at an angle equal to the planet inclination. Hence the 
position angle of the planet orbit can be retrieved. 

We conclude that in most cases, stellar activity will not 
be the dominant factor for the planet detectability in astrome- 
try. The dominant factor will rather be the instrumental preci- 
sion. Note that if present, other larger planets will also induce 
noises, which are not considered in this paper, but the series of 
blind tests conducted to estimate the detection capabilities of 
Earth-like planets in multiple systems by space-borne astrom- 
etry dTraub et al. 20101 1 have showed that this noise can be well 
circumvented. 



3.2. Comparison to RV data 

The situation is then different from the RV data, as the rms of the 
RV due to spots and bright structures only (i.e. convection not 
taken into account) during resp. the entire cycle, the low activity 
period and the high activity periods are 3, 1 and 4 times larger 
than the amplitude of the Earth RV signal (see Table [T}. When 
taking convection into account, the coiTesponding ratio become 
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